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THEORIES OF FLOW SIMILITUDE 

By A. F. Zahm 



SUMMARY 

The laws of comparison of dynamically similar fluid motions are derived by three different 
methods based on the same principle and yielding the same or equivalent formulas. , In this 
report prepared for publication by the National Advisory Committee for Aeronautics, in 
June, 1927, are outlined the three current methods of comparing dynamically similar motions, 
more especially of fluids, initiated respectively by Newton, Stokes (or Hel m holtz), and Rayleigh. 
These three methods, viz., the integral, the differential, and the dimensional, are enough alike 
to be studied profitably together. They will presently be treated in succession then compared. 

INTRODUCTION 

Geometrically similar figures, — If two figures are geometrically similar, they have a constant 
scale ratio 

Ifl^a (1) 

where 1, li are any two homologous lengths. If x, Xi etc., are homologous point coordinates 
for the figures, xfxi = y/yi = z/Zi = a. 

Geometrically similar motions. — Two similar configurations perform geometrically similar 
motions when their homologous points trace similar paths in proportional times; that is, in 
times tj h having any arbitrary ratio 5, the same for all homologous path segments. Thus i>, i\ 
being corresponding path speeds, 

t/h - 1 vlv x = lUflJ = a/5 vth = ViU 1 = l^iW = a[b 2 (2) 

where Ijli is the scale ratio of homologous moving parts, path segments, radii of curvature, etc. 
Since by (2) the ratio liv-[li\ 2 of accelerations normal to the path elements equals bjh along 
them, the resultant accelerations, j } j ± bear the same ratio and are alike directed. The con- 
stant ratios lfl u tft 1} vfv u v/ih all may be different; only two can be independent, as (2) shows. 1 

Dynamically similar systems. — Let the homologous elements of two similar configurations 
in similar motion be masses m, mi having the constant ratio 

m/V/ti = c = p£ 3 /Vi 3 (3) 

p, pi being their densities; then, to keep their motions similar, all corresponding impressed 
forces must be in constant ratio and like direction. 2 For since these elements have resultant 
accelerations j, ji cc t^/Z, Vifl^ their resultant impressed forces R, R t have the ratio 

= mj/m 1 j l = pZV/pA V (4) 

which is constant throughout, since p/p 1? lfl u vfvi are so. Further, the accelerations j f j x are 
alike directed: so then must be R, Ri. So, too, the corresponding forces on large homologous 

i Were the paths similar irrespective of describing time, the motions still would be geometrically similar, but not as de6ned. here and in usual 
writings on similitude. The geometrically similar motions here treated are kinematically similar because they trace similar paths in proportional 
times. 

i That is, their magnitudes are in constant ratio and their lines of action similarly located in the two systems, though the systems themselves 
may be neither simultaneous nor alike oriented in space. 
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parts must be in constant ratio and like direction, as appears on compounding those on their 
constituent elements. 3 Also by the argument for (4) the constituents P, Q, P u Qi/etc, of 
R t R u such as weight, pressure, friction, etc., must be in constant ratio and like direction, viz, 
P to P 1} Q to Q u etc. In fact for homologous elements they are concurrent and have similar 
force polygons. Hence 

BfB l = p1 2 v 2 /pxIM 2 = P/Pi = QfQi - etc _ . (4 t ) 

Such systems are dynamically similar and have (1), (2), (4) as their conditions or criteria of 
similarity. 

By (4) when p/p u ljl u vfv t are assumed constant B/B x is found constant. So, too, if p/pi, 
l/l u R/R 1 are constant, vfvi is constant, and the motions are similar. Fixing either three of these 
ratios determines the fourth. Thus, premised initial similarity, similar mass systems in similar 
motion are similarly forced; conversely similar mobile mass systems similarly forced similarly 
move. In either case the systems are dynamically similar. 

Summation of impressed forces, — The resultant forces B } B x at homologous elements have 
the components 

mj x =F X + <? x '~-fetc.1 
m I ;i x =P lx -hQ 1 r+etc.J 

with like values for the y f z directions. These equations may be compared with (17), where 
the magnitudes of P, Q, etc., not merely their ratios, have definite expression; also with (13), 
where the magnitudes have only proportionate expression. 

DYNAMICALLY SIMILAR FLOWS 

A) NEWTONIAN OR INTEGRAL METHOD 

Definition. — Fluid streams that everywhere satisfy (1), (2), (4) are dynamically similar 
systems, with similar flow fields and boundaries; hence are comparable in their corresponding 
characteristics. 

Classification of chief force ratios. — As before, the ratio of the acceleration forces on homol- 
ogous parts of such systems must be the same throughout and must equal severally the ratios 
of the corresponding impressed forces. The following table exhibits the chief ratios of present 
interest. Their proof follows the table. For all homologous elements the ratios g[g u pfp lT 
p/Hi are assumed constant, p. denoting viscosity. 



TABLE I 

RATIO OF CORRESPONDING FORCES ON HOMOLOGOUS FLUID ELEMENTS 



Ratio of acceleration 
forces mjfmiji 


Ratio of impressed forces 


Gravitational 
mg/migi 


Pres^al, Pg/vgf 


Viscous 


pWtpilW 


gpP/giPih 


fpl 2 v 2 /pih 2 vi 2 , for incompressible fluid 

U^Mk 2 , for elastic fluid.. 


^fdvipikvu 






.. 



Proof of force ratios. — The ratio in column 1 has been proved; that in column 2 is obviously 

true. 

To prove column 3, the pressure force on any small volume of frictionless fluid, being pro- 
portional to volume times along-stream pressure gradient, varies as V.bp/bl, as is well known, 
where bpfblccb(pv 2 )/bl. Hence for p constant the resultant pressure force varies as pPv 2 ; and 
for p variable bpfdl = tcfp.dpfbl, by hydrostatics; that is, the pressure force varies as d 2 , where 
k is the bulk modulus. One recalls that k/k x = pc 2 /piC : 2 where c } c x are the speeds of sound in the 
fluids under the actual working conditions. 



a Newton, reference 1, proves this theorem verbally without using symbols, A different symbolic treatment is given by Sir Richard 
Glazebrook in reference 2. 
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If dv/dl is the rate of distortion in any fluid element, the entailed force on it varies as 
l 2 .fidvfdl oz $v; hence the ratio in column 4. 

Examples of similar How conditions. — Granted kinematic similarity, when the impressed 
forces are as in Table I the general conditions (4i) for dynamic simil arity are 



pZV _ gpl z yip 



pl 2 v 2 /pili% 2 } for incompressible fluids 1 



(6) 



PihW gipii nJiVi 1 k1 2 [k^, for elastic fluids j 

where only the ratios of predominant forces are to be retained. A few examples will illustrate. 

(a) Thus, if weight is the only dominant impressed force, the motions are dynamically 
similar when the first ratio in (6) equals the second, viz, when 

gl/v 2 = g 1 l l fv 1 2 (7) 

which is the well-known Reech and Froude "law of corresponding speeds." 

(j8) If weight and elasticity are negligible, the first ratio in (6) is equated to the third, 
giving 

. (8) 

which is the familiar Reynolds's condition for similarity of motion of fluids. It applies to the 
motion of airships, submarines, skin friction planes, fluids in pipes, etc. 

{'7) If there is considerable compression, while gravity and friction are negligible, the 
first term is equated to the lower fourth, giving 

- = ^_____ ____ ._„..(9) 

which is Booth and Bairstow ? s condition for similarity. 

(5) If g, fx, c all are important, conditions (7), (8), (9) must coexist; if all are negligible, 
(6) gives pl 2 t?fpili% 2 = plVfpili 2 Vi 2 , that is, all flows with similar boundaries are similar, whatever 
the densities and velocities. 4 

Reactions in similar flows. — If P, Pi are corresponding reactions of a craft and its model 
under conditions (7), P/P l =gpl z fg lPl l 1 z , whence 

P=N igp lK- — - (10) 

where Ni = P\lgipJ<i is a dimensionless coefficient, say, given by model tests. 

If g, k are negligible, conditions (S) obtain, and P/Pi = pLvjiL^ x v u or P/Pi = plVfpiliV 2 , whence 

P=N 2 !dv, or P = OplV (11) 

where N2 = Pi/f*ilii\ } 0=PiIp\1 2 v 2 } both dimensionless coefficients. 

If g } p are negligible, and compression important, P/P^^V^i 2 ; hence 

P=2V 3 /cL__...__ (12) 

where i\ 7 3 = Pi/^Z 1 2 , and conditions (9) prevail. 

Let the Ps be all lifts or all drags or other like directed forces. Then, if g, fx, k all are impor- 
tant together, the total of such reactions on the craft is 

R = N igp l* + N 2 ult>+ Nspc 2 P= plWfiglfv 2 , v/lv } c/»)=AW-- - 5 (13) 

got by summing (10), (11), (12), using *:=pe 2 , then factoring off pZ¥\ One notes that (13) can 
be written: Total reaction = gravitational ^frictional-f pressural. 

The validity of (13) was premised on dynamic similarity of motion of the craft and its 
model, as defined by the simultaneous conditions (7), (8), (9). That is, 

N ig l/v 2 ^ N 2 v/lv^ N z c 2 /v^f(gl/i^ vflv, c/v) =f($ 1 l l [v 1 2 , vjl&u c 1 fv l ) = N. . (14) 

* (7) has the alternative form g!i=Qi}j\: (S) the alternative ffp=filpi, where , fi, p t pi are corresponding frictions and pressures per unit area of 
homologous surface elements 8S,5Sl With fji L top, p t the resultant stresses have slopes ftp, f if pi to the normals at 5 S, 5 Si. The "laws" (7), (8) 
(9) are but corollaries of (6) or (40 . 

* More conventionally one writes R=plWfi ts*fgl, rfto, sic). 
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Alternatively (13) can be written 

B = iilvf(gl/v 2 } v/lv, c/v) = N'txh- 



-(130 



found by factoring off fdv and rearranging the results 

Writers sometimes say (13) shows that R varies as pl 2 v 2 ; they can as well say Roc^lv by (13!). 
The first statement is true for conditions making/ constant; the second for/ 7 constant. (Fig. 1 .) 

Dynamic scale and scale effect— II N x ( = Pifgml*), found from a model test ; is plotted against 
gJi/vi 2 i^gl/v 2 ) the graph is directly applicable to computing the full-scale reaction (10); simi- 
larly for the graph of A r 2 against v/lv and N z against c/v. 

In such plots the dimensionless argument, say, v/lv, is treated as a single independent 
variable. The graph is the same whether v varies alone or I alone or v alone, or if two or three 
vary together. If N 2 varies as (lv/v) n } it varies as l n , v n , v~ n . The effect, for instance, of varying 
I can be learned by varying v or v in the model test, and so for JVi, N 3 . 

One calls the independent dimensionless argument gl/v 2 the dynamic scale for the motion 
(a), and the variation of N x with scale the scale effect. Similarly v/lv is the dynamic scale for 




!0' s 



to-* t 2 5 (0 .10* 

Reynolds Number, Vdl» 



F IQ . i— Drag coefficients Cd^DIpVW, Cd'=DIiiV<1, plotted against Vdfr for a sphere in uniform translation 
through a viscous fiuid. Data given in N\ A. C. A. Report No. 253. Either graph can he plotted from the 
other, since Cs>'iCD=Vd{v 

motion (/?), and l/c for (7). No doubt the term "scale effect" originally meant the effect of 
changing the linear s.cale ratio lfl u then was extendedto mean the effect of changing some more 
complex argument, such as gl/tP, v/lv, c/v, etc., now "called the dynamic scale. The simpler 
scale ratios l/l u t/t u p/p u etc., are called scales of length, time, pressure, etc. 

The more complex reaction (13) is a function of three dynamic scales, shown in parentheses. 
The scale effect here is the variation of iV or R with one or more of the dynamic scales, or inde- 
pendent arguments gl/v 2 } v/lv } c/v. But for the particular case g f v } c = o, as for a perfect liquid 
unaffected by gravity, N h N 2} iV 3 are constant and have straight-line graphs when plotted 
against their scales. Then (13) gives R = const, times pl 2 v 2 . 

Generally, therefore, for dynamically similar fluid motions a dynamic scale is any one of the 
independent dimensionless arguments in the formula for the fluid reaction; the scale effect- is 
the variation of such reaction due to variation of the arguments. 

Arbitrary and derived scale ratios. — As seen in the.introduetion, geometric similarity requires 
one constant scale ratio, say, lfl h for length; kinematic similarity two scales, say, lfl u tft u for 
length and time; dynamic similarity three, say, lfl u t/t u p/p x . From these many others may be 
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derived, say, vfv u vfvi, pfpi, vjv u etc., in case of fluids. Any three can be taken as determinative, 
then combined to form derived ratios, as exemplified in (2). For dynamic similarity we arbi- 
trarily chose, at the outset, ljl u tjt u m\m u because length, time, and mass usually appear as 
basic. For the same service with fluid systems Helmholtz (reference 3) takes p/p u vjv^ vfv u 
while other writers choose still other scales as fundamental. 

Thus for geometrically similar fluid motions Helmholtz, assuming 

Pifp=r • = g uJu^Vi/v=w 1 fw=^n (15) 

as given constants, thence derives the further ratios 

Xifx = yify = Zifz = q/n f x /f = g_[n 2 p^p = n 2 rp + const (16) 

for use in comparing the differential equations of motion of the two fluids. 

(B) DIFFERENTIAL METHOD 

The conditions (6) for dynamically similar motion of two fluids can also be derived from the 
standard differential equations of motion of such fluids, viz, from 



-dp «, 1 bd 



-(17) 



with like expressions for the y, z and y u z 1 directions. 6 For if the motions are dynamically 
similar corresponding terms, all being forces, must have the same ratio, 

pit, pg x c 2 bp/dx t±A 2 u _ ydd/dx 
Piu~ Pi^ix - ^ 2 dpi/d^" ^A 2 ^ mbdtldxi 

Expressing these ratios in finite dimensions gives 

pXjU 2 pg x _ px x c 2 _ px-ihi 
PiXUi 2 ~ Pigi x ~ PiXd 2 ~ ixix?u x 

which multiplied by 2?fxi z become the relations (6) for the x direction, viz, 

px 2 u 2 _ gpx 2 _ px 2 c 2 _ pxu 
Pi^iW ~giPiXi "pi^iV ~mx&i 

Thus the differential method yields the same result as the Newtonian. It is Newton's 
method in Stokes's shorthand, except that Stokes would first write the forces, then their ratio; 
Newton would write their ratio directly. But to write their ratio one must know approximately 
their nature and analytic expression. 

Helmholtz reverses the above argument. Assuming the relations (15), (16), he says they 
transform the first of (17) into the second, omitting the g terms. Hence he infers that model 
data serve to predict the hydrodynamic behavior of full-scale craft when the relations (15), 
(16) are maintained. 

Diverse and sundry treatments of this topic are found in references 8, 9, 10. 

(C) DIMENSIONAL METHOD 

Mechanical units. — In mechanics three measuring units, say, of length, mass, time— L 3 
If, T, arbitrarily taken as fundamental, are combined in various powers to form other kinds 
called derived units, such as U=AL*H Y T % i A being constant, and x, y, z > 1. Table II 
illustrates. These two classes of units, viz, fundamental and derived, serve to measure 
mechanical quantities of every kind, such as length, speed, torque, etc. Thus any mechanical 
quantity f?, if a function of n others, all differing in kind, can be written 

■R = 2MQ l *Q 2 b —-Q™-„ (18) 

where if, a t b m are pure numbers, and Q may involve either fundamental or derived units. 

Table II gives examples. 

+ ff-^+^+^t etc., for Aim,*. 
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Basic formula. — To be valid for all unit systems (18) must be dimensionally homogeneous. 7 
Then it can be written, if LR] = [£i a <>2 b — ~£n m ] = LP], 

B-2tfP-„- - - (19) 

where iVis a pure number. This is the basic formula of dimensional theory, and is most general 
when the Ps are all the separate independent Q products that can be formed having the dimen- 
sions of R. B 

Homogeneous products. — To find these P products/vve first multiply any Q triad, say Q x Q y Q\ 
by each remaining Q successively, and equate to [R] the dimensions of each resulting product. 
The ensuing n-3 equations, with tentative exponents a } 5, c, are 

[Pi] = {Qi & Q2 h Qz c Q^\R] [P2}-[Qi 6 Q2 e Q^Q,} = [R], etc.™- (20) 

Now replacing [R] and [Q] by their values in L } M } T and solving (20) for a, b } c— gives P u 
P 2j The following example illustrates. For a rigorous analysis special works on dimen- 
sional theory may be consulted (references 4, 5). 

Reactions in similar flows. — If the reaction R of. a body in a fluid stream depends solely 
on p, l t v, g, c } or density, size, speed, weight, viscosity, elasticity, all the separate independent 
products having the dimensions of R possible to make with them amount to 6-3, say, P u P 2 ? P». 
To form these, we take any triad p x l y v z of the six independent quantities and multiply it succes- 
sively by the remaining ones g, jx } c } giving 

P, - p*T°v\g P 2 =-p W.ju P 5 = p"ZV.c 

and equate the dimensions of each product to [K\ = [ML/T 2 .] The first yields, by Table II, 

(M/L*)\L h .(L/Ty.L/T 2 = ML/T 2 

each unit having the same aggregate exponent in both terms. On equating the indices of 
L, M y T successively this gives 

-3a+6 + c+l=l a = l_ -c-2=-2 

Thus a= 1, b = 3, c = 0, whence P u = pVg. A like procedure gives d=*0, 6=1 =/, whence P 2 = lv]x. 
Similarly g=l } & = 2, whence P z = pl 2 vc. 
By (19) the reaction now is 

R~N X pl z g + N 2 iilv+ N zP l 2 vc = f>Wf{gllv\ v/lv, c/v) . .-(13) 

which is a general resistance equation for the specified dynamical conditions, viz, that R is a 
function solely of p, I, v, g } fx, c. From this (130 also is found as before. 

By (13), if the arguments in parentheses are given any specific values the same for model 
and full scale, / is the same for both; hence 

R^NplV 

where N=f(g x lifv 2 y vifl x v u the same as by Newton's method. 

COMPARISON OF THE THREE METHODS 

In the foregoing text the same criteria for dynamical similarity in two flow systems were 
found by three different methods of analysis — the Newtonian, the differential, and the dimen- 
sional. In each the physical quantities governing the flow were premised from experience. 
Thence were found the ratios of corresponding impressed forces of each kind on homologous parts 
of the fluid, viz, weight ratio, pressure ratio, etc. These ratios, by definition of dynamic simi- 
larity, must each equal the ratio of the resultant acceleration forces on those parts; viz, the ratio 

By Newton's method we directly equated the ratio of these acceleration forces to the 
several ratios of corresponding impressed forces, thus obtaining specific conditions for dynamical 



" That is, all terms of (18) must comprise the same fundamental units, each having a constant aggregate exponent throughout the equation. 
f Dividing (19) by R gives $ On, tt 2 v\) =0, the being dimensionless products. This is Buckingham's * theorem (reference 4). 



THEORIES OF FLOW SIMILITUDE 



193 



similarity, and formulas for the reactions of any fluid system in terms of those of its modeL By 
the second method we first wrote the differential equations for the two fluid motions assumed 
dynamically similar, then equated the ratios of corresponding terms, thus obtaining the same 
result as before. By the third method we first equated the unknown reaction B on one fluid 
element to the sum of all the terms we could form from the flow-governing quantities arranged in 
power products each having the dimensions of R. Doing the same for the homologous element, 
then taking the ratio of the R forces, gave the same reaction formula as found by the other methods. 

At first sight the dimensional process seems to be a routine algebraic operation requiring 
less knowledge than is needed for the two other methods. In reality all three demand adequate 
judgment of the kind of physical quantities governing the motion, and their comparative 
importance. In all three cases the assumed physical agencies are the same, the terms in the 
dynamic equations are analogous, and the final working formulas are the same or equivalent. 
In all, the derived working formula contains a dimensionless coefficient that is not deduced 
theoretically, but is to be found from model tests, then applied to full-scale apparatus operating 
under dynamically similar conditions. In all, the "laws of comparison" are merely expressions 
of equality of like dynamic scales, viz, equality of the ratio of the acceleration forces to the 
corresponding ratios of the dominant impressed forces. 

TABLE II 

QUANTITIES EXPRESSED IN BASIC UNITS OF LENGTH, TIME, AND MASS, L, T, M 

The "dimensions" of a physical quantity are the degrees of the fundamental units in its formula. Thus the dimensions of an acceleration, 
which are symbolized by [L T~ z ] f are 1 in length and —2 in time. Commonly the brackets are omitted from such simple L } T f M expressions 
not containing other'symbols. The dimensions of a force are MLT^- I viz, 1 in mass, 1 in length, —2 in time; the dimension of an angle, a sine, 
cosine, tangent, etc., is L,L~ L , that is zero. Logarithms in physical equations operate only on dimensionless quantities, such as pure numbers 
or ratios of like physical quantities; hence are dimensionless. 

A derived unit, being formed of powers of fundamental units, has the form V= AL*M*T* t with dimensions j>J_V2X Thus a force F=ms[fl= 
A ML T~-, where m, s, t are mass, length, time in any convenient units. Its dimensions are written [F]= [mslt-\=lA ML T~ l \ — ML T~*. 

In homogeneous equations all terms have the same dimensions, that is, the same aggregate exponent for each basic unit. Thus in the last equa- 
tion of Table II each term has the dimensions ML* T~'-: for [pv*) = ML^.L 1 T<= ML* T~\ and UiB ( 'a] = ML* T*.LT*IL= ML* T*. In the 
familiar projectile formula c=*gt-\-c,\g([=LT-* r T*=LT-i=lvl=[c], where c (=e e ) is a velocity. 



Kind of quantity 


Symbol. 


Formula 


Dimensions of each term 




Derived units. 







Area, surface 

Volume 

Angle 

Linear velocity _ 

Linear acceleration 

Angular velocity 

Angular acceleration 

Density 

Force, thrust 

Torque, moment 

Pressure, friction 

Work, energy, potential. 

Power, activity 

Viscosity 

Kinematic viscosity 

Flux of fluid 

Velocity potential 



S=P _ 

r=Z* 

e = 5 /r = arc -f- radius _ _ . 
u = x[t= dxf dt=x 

. j=ult=duldt=d-zldP. 
u = 8/i= d6jdi= du/d#_ _ 

a=c*it=du[dt=e[t 2 

p = 7n/r=mfl z 

F=mj=ms[P 

Q-Ft 

p=F/S s f=F/S 

W=Fs 

P=Fu 

^=J*-du/dy=f/o} 

y=Mf'p 

6=fq a dS 

<p= —fqtds 



LK 

LK 

L°. 

LT~K 

LT-*. 

T~ l . 

r p-i 

ML" 3 . 

MTr i T - K 
MI?T~*. 

Mlr l T~ l . 

L % T~ X . 



Geometrical and mechanical equations, R — 'ZNP 



Length of catenary 

Area of ellipse 

Volume of frustum of cone- 



Period of simple pendulum 

Mutual attrac. of two particles, _. 
Strength of line source _. 



<? for source-sink in plane stream. . 

Acceleration of viscous particle. 
Nose pressure on falling droplet 



= £ ( e x/a_ e -x/c) = c sin A ^) — 



S = 7cab- 
r=-|xA(r 2 H-rr'-r-r /2 )_ 



t=2^Ug _ 

F=Kmm l fr 1 , where [kI^U/MT 2 .. _ 
m — 2-aaq n 

<p=iix+(r log^r, where [<r] = I?/T- 

u~ — dp/p bx-rvA 2 u 

p tt =pi' 2 /2-f ISfie/a 



= L.LK 



= [L(I? + LL+L*)]. 



&T-\ =L.LT~K 

L*T~\ =[LT T - I .X-FL 2 7 T - 1 log L% 
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SYMBOLS USED IN TEXT 



?, i Homologous lengths in similar figures. 

x, Xi, etc Homologous coordinates. 

t, fr— _ Times .of tracing homologous paths. 

v } vi- Corresponding path velocities. 

v v\ Corresponding path accelerations. 

j 9 ji Corresponding total accelerations. 

m, mi Homologous masses. . ... 

a, b, c Arbitrary numerical ratios l/h, tjti, ml mi) also tentative exponents. 

p t ju, v, k Density, viscosity, kinematic viscosity, bulk modulus; ditto for p u n u v u k x . 

c=V^/p Speed of sound in elastic fluid; ditto for ci. 

R, R t Resultant forces on mass elements m, mi. 

P, Q f etc Components of R; P u Gu etc., components of R\. 

p f f Pressure and friction, per unit area; ditto for p u fa. 
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